Evidence for oscillation of atmospheric neutrinos 
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We present an analysis of atmospheric neutrino data from a 33.0 kiloton-year (535-day) 
exposure of the Super-Kamiokande detector. The data exhibit a zenith angle dependent 
deficit of muon neutrinos which is inconsistent with expectations based on calculations of 
the atmospheric neutrino flux. Experimental biases and uncertainties in the prediction 
of neutrino fluxes and cross sections are unable to explain our observation. The data 
are consistent, however, with two-flavor v ^ <-> v T oscillations with sin 2 28 > 0.82 and 
5 x 10~ 4 < Am 2 < 6 x 10~ 3 eV 2 at 90% confidence level. 
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Atmospheric neutrinos are produced as decay 
products in hadronic showers resulting from col- 
lisions of cosmic rays with nuclei in the upper at- 
mosphere. Production of electron and muon neu- 
trinos is dominated by the processes ir + — > /i + + 
followed by /i + — > e + + + v c (and their charge 
conjugates) giving an expected ratio (= v^/vc) of 
the flux of + to the flux of v c + v c of about 
two. The Vfi/fe ratio has been calculated in detail 
with an uncertainty of less than 5% over a broad 
range of energies from 0.1 GeV to 10 GeV Jl|,||. 

The v^/ve flux ratio is measured in deep under- 
ground experiments by observing hnal-state lep- 
tons produced via charged-current interactions of 
neutrinos on nuclei, v + N — > I + X. The flavor of 
the final state lepton is used to identify the flavor 
of the incoming neutrino. 

The measurements are reported as R = 
e) data/ (/•*/ g) mc i where fi and e are the num- 
ber of muon-like (/x-like) and electron-like (e-like) 
events observed in the detector for both data and 
Monte Carlo simulation. This ratio largely cancels 
experimental and theoretical uncertainties, espe- 
cially the uncertainty in the absolute flux. R = 1 is 
expected if the physics in the Monte Carlo simula- 
tion accurately models the data. Measurements of 
significantly small values of R have been reported 
by the deep underground water Cherenkov detec- 
tors Kamiokande || : 4], 1MB ||, and recently by 
Super-Kamiokande [fiQ. Although measurements 
of R by early iron-calorimeter experiments Frejus 
H and NUSEX [|| with smaller data samples were 
consistent with expectations, the Soudan-2 iron- 
calorimeter experiment has reported observation 
of a small value of R 10 1 . 

Neutrino oscillations have been suggested to ex- 
plain measurements of small values of R. For a 
two-neutrino oscillation hypothesis, the probabil- 
ity for a neutrino produced in flavor state a to be 
observed in flavor state b after traveling a distance 
L through a vacuum is: 



Pn 



sin 29 sin 



T.27Am 2 (eV 2 )Z(km) 
. E v (GcV) 
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where E v is the neutrino energy, 9 is the mixing 
angle between the flavor eigenstates and the mass 
eigenstates, and Am 2 is the mass-squared differ- 
ence of the neutrino mass eigenstates. For detec- 
tors near the surface of the Earth, the neutrino 
flight distance, and thus the oscillation probabil- 
ity, is a function of the zenith angle of the neu- 
trino direction. Vertically downward-going neutri- 



nos travel about 15 km while vertically upward- 
going neutrinos travel about 13,000 km before in- 
teracting in the detector. The broad energy spec- 
trum and this range of neutrino flight distances 
makes measurements of atmospheric neutrinos sen- 
sitive to neutrino oscillations with Am 2 down to 
10 -4 eV 2 . The zenith angle dependence of R mea- 
sured by the Kamiokande experiment at high en- 
ergies has been cited as evidence for neutrino os- 
cillations [0. 

We present our analysis of 33.0 kiloton-years 
(535 days) of atmospheric neutrino data from 
Super-Kamiokande. In addition to measurements 
of small values of R both above and below ~1 GeV, 
we observed a significant zenith angle dependent 
deficit of /i-like events. While no combination of 
known uncertainties in the experimental measure- 
ment or predictions of atmospheric neutrino fluxes 
is able to explain our data, a two-neutrino oscil- 
lation model of «-> v x , where v x may be v r or 
a new, non-interacting "sterile" neutrino, is con- 
sistent with the observed flavor ratios and zenith 
angle distributions over the entire energy region. 

Super-Kamiokande is a 50 kiloton water Cheren- 
kov detector instrumented with 11,146 photomul- 
tiplier tubes (PMTs) facing an inner 22.5 kiloton 
fiducial volume of ultra-pure water. Interaction 
kinematics are reconstructed using the time and 
charge of each PMT signal. The inner volume is 
surrounded by a ^2 meter thick outer detector in- 
strumented with 1885 outward-facing PMTs. The 
outer detector is used to veto entering particles and 
to tag exiting tracks. 

Super-Kamiokande has collected a total of 4353 
fully-contained (FC) events and 301 partially- 
contained (PC) events in a 33.0 kiloton-year ex- 
posure. FC events deposit all of their Cherenkov 
light in the inner detector while PC events have 
exiting tracks which deposit some Cherenkov light 
in the outer detector. For this analysis, the neu- 
trino interaction vertex was required to have been 
reconstructed within the 22.5 kiloton fiducial vol- 
ume, defined to be > 2 m from the PMT wall. 

FC events were separated into those with a sin- 
gle visible Cherenkov ring and those with mul- 
tiple Cherenkov rings. For the analysis of FC 
events, only single-ring events were used. Single- 
ring events were identified as e-like or /^-like based 
on a likelihood analysis of light detected around 
the Cherenkov cone. The FC events were sep- 
arated into "sub-GeV" (E vis < 1330 MeV) and 
"multi-GeV" {E ms > 1330 MeV) samples, where 



2 



E V i S is defined to be the energy of an electron that 
would produce the observed amount of Cherenkov 
light. E V i a = 1330 MeV corresponds to ~ 1400 
MeV/c. 

In a full-detector Monte Carlo simulation, 88% 
(96%) of the sub-GeV e-like (//-like) events were v e 
(i/fj.) charged-current interactions and 84% (99%) 
of the multi-GeV e-like (/i-like) events were v c (v^) 
charged-current interactions. PC events were esti- 
mated to be 98% charged-current interactions; 
hence, all PC events were classified as /x-like, and 
no single-ring requirement was made. Table | sum- 
marizes the number of observed events for both 
data and Monte Carlo as well as the R values for 
the sub-GeV and multi-GeV samples. Further de- 
tails of the detector, data selection and event re- 
construction used in this analysis are given else- 
where fpPl. 



Data Monte Carlo 


sub-GeV 








single-ring 


2389 


2622.6 




e-like 


1231 


1049.1 




/i-like 


1158 


1573.6 




multi-ring 


911 


980.7 




total 


3300 


3603.3 




R = 


0.63 ± 0.03 (stat.) 


± 0.05 (sys.) 


multi-GeV 








single-ring 


520 


531.7 




e-like 


290 


236.0 




//-like 


230 


295.7 




multi-ring 


533 


560.1 




total 


1053 


1091.8 




partially-contained 


301 


371.6 




Rfc+pc = 


0.65 ± 0.05 {stat.) 


± 0.08 (sys.) 



TABLE I. Summary of the sub-GeV, multi-GeV 
and PC event samples compared with the Monte Carlo 
prediction based on the neutrino flux calculation of 
Ref. [2]. 



We have measured significantly small values of 
R in both the sub-GeV and multi-GeV samples. 
Several sources of systematic uncertainties in these 
measurements have been considered. Cosmic ray 
induced interactions in the rock surrounding the 
detector have been suggested as a source of e-like 
contamination from neutrons, which could pro- 
duce small R values pj| , but these backgrounds 
have been shown to be insignificant for large water 
Cherenkov detectors [^2|. In particular, Super- 
Kamiokande has 4.7 meters of water surrounding 
the fiducial volume; this distance corresponds to 
roughly 5 hadronic interaction lengths and 13 ra- 
diation lengths. Distributions of event vertices ex- 
hibit no excess of e-like events close to the fiducial 
boundary (|J^]. 

The prediction of the ratio of the flux to the 
v e flux is dominated by the well-understood de- 
cay chain of mesons and contributes less than 5% 
to the uncertainty in R. Different neutrino flux 
models vary by about ±20% in the prediction of 
absolute rates, but the ratio is robust |n|. Uncer- 
tainties in R due to a difference in cross sections 
for v c and have been studied ]L4| ; however, lep- 
ton universality prevents any significant difference 
in cross-sections at energies much above the muon 
mass and thus errors in cross-sections could not 
produce a small value of R in the multi-GeV en- 
ergy range. Particle identification was estimated to 
be <: 98% efficient for both ^t-like and e-like events 
based on Monte Carlo studies. Particle identifi- 
cation was also tested in Super-Kamiokande on 
Michel electrons and stopping cosmic-ray muons 
and the /i-like and e-like events used in this analy- 
sis are clearly separated §. The particle identifi- 
cation programs in use have also been tested using 
beams of electrons and muons incident on a water 
Cherenkov detector at KEK |lf|. The data have 
been analyzed independently by two groups, mak- 
ing the possibility of significant biases in data se- 
lection or event reconstruction algorithms remote 
Other explanations for the small value of R, 
such as contributions from nucleon decays 16 1, can 
be discounted as they would not contribute to the 
zenith angle effects described below. 

We estimate the probability that the observed 
/i/e ratios could be due to statistical fluctuation is 
less than 0.001% for sub-GeV R and less than 1% 
for multi-GeV R. 

The //-like data exhibit a strong asymmetry in 
zenith angle (6) while no significant asymmetry 
is observed in the e-like data M. The asymme- 
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try is denned as A = (U - D)/{U + D) where 
U is the number of upward-going events (—1 < 
cos 9 < —0.2) and D is the number of downward- 
going events (0.2 < cosG < 1). The asymmetry is 
expected to be near zero independent of flux model 
for E v > 1 GeV, above which effects due to the 
Earth's magnetic field on cosmic rays are small. 
Based on a comparison of results from our full 
Monte Carlo simulation using different flux mod- 
els Jl]]^] as inputs, treatment of geomagnetic effects 
results in an uncertainty of roughly ±0.02 in the 
expected asymmetry of e-like and /i-like sub-GeV 
events and less than ±0.01 for multi-GeV events. 
Studies of decay electrons from stopping muons 
show at most a ±0.6% up-down difference in Che- 
renkov light detection fl7|| . 

Figure [j] shows A as a function of momentum for 
both e-like and /it-like events. In the present data, 
the asymmetry as a function of momentum for e- 
like events is consistent with expectations, while 
the /i-like asymmetry at low momentum is consis- 
tent with zero but significantly deviates from ex- 
pectation at higher momentum. The average angle 
between the final state lepton direction and the in- 
coming neutrino direction is 55° at p — 400 MeV/c 
and 20° at 1.5 GeV/c. At the lower momenta in 
Fig. |l|, the possible asymmetry of the neutrino flux 
is largely washed out. We have found no detector 
bias differentiating e-like and /i-like events that 
could explain an asymmetry in //-like events but 
not in e-like events [0. 

Considering multi-GeV (FC±PC) muons alone, 
the measured asymmetry, A = —0.296 ± 0.048 ± 
0.01 deviates from zero by more than 6 standard 
deviations. 

We have examined the hypotheses of two-flavor 
<-> v c and «-> v T oscillation models using a 
X 2 comparison of data and Monte Carlo, allow- 
ing all important Monte Carlo parameters to vary 
weighted by their expected uncertainties. 

The data were binned by particle type, momen- 
tum, and cos 9. A \ 2 is defined as: 



cosQ,p 



DATA 



MC, 



v* 2 (2) 



where the sum is over five bins equally spaced 
in cos 9 and seven momentum bins for both e-like 
events and /i-like plus PC events (70 bins total). 
The statistical error, er, accounts for both data 
statistics and the weighted Monte Carlo statistics. 
Ndata is the measured number of events in each 



bin. Nmc is the weighted sum of Monte Carlo 
events: 



Lb at a 
N mc = — x 



MC 



E 



(3) 



MC events 



£-data and Cmc are the data and Monte Carlo 
live-times. For each Monte Carlo event, the weight 
w is given by: 



x < 



(l + a)(££/£ ) 5 (l4 

/ e>M (sin 2 20,Am 2 ,(l 

(l-ft/2) 
(1+&/2) 

(1 - Pm/2) 



7/ Sim cos9) x 

+ X)L/Eu) 

sub-GeV e-like 
sub- GeV /x-like 
multi-GeV e-like 



(1 + A„/2)(l - § £jp) multi-GeV /i-like 



(1 + /W2)(1 



E ) 



PC 

(4) 



E l v is the average neutrino energy in the i th mo- 
mentum bin; E is an arbitrary reference energy 
(taken to be 2 GeV); rj s (r) m ) is the up-down un- 
certainty of the event rate in the sub-GeV (multi- 
GeV) energy range; Npc is the total number of 
Monte Carlo PC events; iV M is the total number 
of Monte Carlo FC multi-GeV muons. The fac- 
tor f e ^ weights an event accounting for the initial 
neutrino fluxes (in the case of v ^ <-> v e ), oscilla- 
tion parameters and L/E v . The meaning of the 
Monte Carlo fit parameters, a and tj =(f3 s , flm, 
5, p, A, rj s , 7] m ) and their assigned uncertainties, 
<jj , are summarized in Table [n| The over-all nor- 
malization, a, was allowed to vary freely. The un- 
certainty in the Monte Carlo LjE v ratio (A) was 
conservatively estimated based on the uncertainty 
in absolute energy scale, uncertainty in neutrino- 
lepton angular and energy correlations, and the 
uncertainty in production height. The oscillation 
simulations used profiles of neutrino production 
heights calculated in Ref. |l8[| , which account for 
the competing factors of production, propagation, 
and decay of muons and mesons through the atmo- 
sphere. For Vn <-> v e , effects of matter on neutrino 
propagation through the Earth were included fol- 
lowing Ref. [ pOpl] ]. Due to the small number of 
events expected from r-production, the effects of 
t appearance and decay were neglected in simula- 
tions of Vfj, <-» v T . A global scan was made on a 
(sin 2 29, log Am 2 ) grid minimizing \ 2 with respect 
to a, j3 s , j3 m , 5, p, A, rjs and r\ m at each point. 



The best-fit to v„ 



v T oscillations, x 



2 

mm 



65.2/67 DGF, was obtained at (snT 26 
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1.0, Am 2 = 2.2 x 1(T 3 cV 2 ) inside the physi- 
cal region (0 < sin 2 29 < 1). The best-fit val- 
ues of the Monte Carlo parameters (summarized 
in Table |j) were all within their expected er- 
rors. The global minimum occurred slightly out- 
side the physical region at (sin 2 29 = 1.05, Am 2 = 
2.2 x 10~ 3 eV 2 ,xLn = 64.8/67 DOF). The con- 
tours of the 68%, 90% and 99% confidence inter- 
vals are located at Xmin + 2-6, 5.0, and 9.6 based 
on the minimum inside the physical region [jlSfj . 
These contours are shown in Fig. |^. The region 
near \ 2 minimum is rather flat and has many local 
minima so that inside the 68% interval the best-fit 
Am 2 is not well constrained. Outside the 99% al- 
lowed region the x 2 increases rapidly. We obtained 
X 2 = 135/69 DOF, when calculated at sin 2 29 = 0, 
Am 2 = (i.e. assuming no oscillations). 

For the test of «-> v c oscillations, we ob- 
tained a relatively poor fit; Xmin — 87.8/67 DOF, 
at (sin 2 26> = 0.93, Am 2 = 3.2 x 10~ 3 cV 2 ). The ex- 
pected asymmetry of the multi-GeV e-like events 
for the best-fit «-> v c oscillation hypothesis, 
A = 0.205, differs from the measured asymmetry, 
A = -0.036 ± 0.067 ± 0.02, by 3.4 standard devia- 
tions. We conclude that the <-> v e hypothesis is 



not favored. 


Monte Carlo Fit Parameters 


Best Fit Uncertainty 


a overall normalization 


15.8% 


(*) 




8 E v spectral index 


0.006 


(TS = 


0.05 


/3 S sub-GeV fj,/e ratio 


-6.3% 


Cs = 


8% 


fi m multi-GeV \xje ratio 


-11.8% 


Cm = 


= 12% 


p relative norm, of PC to FC 


-1.8% 


CTp = 


8% 


A L/E v 


3.1% 


(J\ = 


15% 


r/ s sub-GeV up-down 


2.4% 


s 

(T n = 


2.4% 


rjm multi-GeV up-down 


-0.09% 


°r, = 


= 2.7% 



TABLE II. Summary of Monte Carlo fit parameters. 
Best-fit values for <-» i/ r (Am 2 = 2.2 x 10~ 3 eV 2 , 
sin 2 26 = 1.0) and estimated uncertainties are given. 
'■"'The over-all normalization (a) was estimated to 
have a 25% uncertainty but was fitted as a free pa- 
rameter. 



The zenith angle distributions for the FC and 
PC samples are shown in Fig. |[ The data are 
compared to the Monte Carlo expectation (no os- 
cillations, hatched region) and the best-fit expec- 
tation for «-> v r oscillations (bold line). 

We also estimated the oscillation parameters 
considering the R measurement and the zenith an- 
gle shape separately. The 90% confidence level 
allowed regions for each case overlapped at 1 x 
10~ 3 < Am 2 < 4 x 10~ 3 eV 2 for sin 2 29 = 1. 

As a cross-check of the above analyses, we have 
reconstructed the best estimate of the ratio LjE v 
for each event. The neutrino energy is estimated 
by applying a correction to the final state lepton 
momentum. Typically, final state leptons with 
p ~ 100 MeV/c carry 65% of the incoming neu- 
trino energy increasing to ~85% at p = 1 GeV/c. 
The neutrino flight distance L is estimated follow- 
ing Ref. fl8|| using the estimated neutrino energy 
and the reconstructed lepton direction and flavor. 
Figure |^ shows the ratio of FC data to Monte Carlo 
for e-like and /i-like events with p > 400 MeV/c as 
a function of L/ E Vl compared to the expectation 
for <-> v T oscillations with our best-fit parame- 
ters. The e-like data show no significant variation 
in L/E v , while the /Li-like events show a significant 
deficit at large L/E v . At large L/E v , the have 
presumably undergone numerous oscillations and 
have averaged out to roughly half the initial rate. 

The asymmetry A of the e-like events in the 
present data is consistent with expectations with- 
out neutrino oscillations and two-flavor v c «-> 
oscillations are not favored. This is in agreement 
with recent results from the CHOOZ experiment 
(22). The LSND experiment has reported the ap- 
pearance of v e in a beam of produced by stopped 
pions [^H . The LSND results do not contradict the 
present results if they are observing small mixing 
angles. With the best-fit parameters for <-> in- 
osculations, we expect a total of only 15-20 events 
from v r charged-current interactions in the data 
sample. Using the current sample, oscillations be- 
tween and v r are indistinguishable from oscil- 
lations between and a non-interacting "sterile" 
neutrino. 

Figure shows the Super-Kamiokande results 
overlaid with the allowed region obtained by 
the Kamiokande experiment The Super- 

Kamiokande region favors lower values of Am 2 
than allowed by the Kamiokande experiment; how- 
ever the 90% contours from both experiments have 
a region of overlap. Preliminary studies of upward- 
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going stopping and through-going muons in Super- 
Kamiokande |f 24 f give allowed regions consistent 



[7] Super-Kamiokande Collaboration, Y. Fukud a, et 
al, submitt ed to Phys. Lett. B (May 1998), |hep 



with the FC and PC event analysis reported in 



t/9805006 



this paper. 

Both the zenith angle distribution of /i-like 
events and the value of R observed in this experi- 
ment significantly differ from the best predictions 
in the absence of neutrino oscillations. While un- 
certainties in the flux prediction, cross sections, 
and experimental biases are ruled out as explana- 
tions of the observations, the present data are in 
good agreement with two-flavor «-> v T oscilla- 
tions with sin 2 29 > 0.82 and 5 x f (T 4 < Am 2 < 
6 x f 0~ 3 eV 2 at 90% confidence level. We conclude 
that the present data give evidence for neutrino os- 
cillations. 
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FIG. 1. The (U - D)/(U + D) asymmetry as a 
function of momentum for FC e-like and /i-like events 
and PC events. While it is not possible to assign 
a momentum to a PC event, the PC sample is esti- 
mated to have a mean neutrino energy of 15 GeV. The 
Monte Carlo expectation without neutrino oscillations 
is shown in the hatched region with statistical and sys- 
tematic errors added in quadrature. The dashed line 
for pi-like is the expectation for Vp. <-> v T oscillations 
with (sin 2 26 = 1.0, Am 2 = 2.2 x 1(T 3 eV 2 ). 
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FIG. 2. The 68%, 90% and 99% confidence inter- 
vals are shown for sin 2 28 and Am 2 for <-> v T 
two-neutrino oscillations based on 33.0 kiloton-years 
of Super-Kamiokande data. The 90% confidence in- 
terval obtained by the Kamiokande experiment is also 
shown. 
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FIG. 3. Zenith angle distributions of /i-like and e-like events for sub-GeV and multi-GeV data sets. Up- 
ward-going particles have cos 6 < and downward-going particles have cos 6 > 0. Sub-GeV data are shown 
separately for p < 400 MeV/c and p > 400 MeV/c. Multi-GeV e-like distributions are shown for p < 2.5 GeV/c 
and p > 2.5 GeV/c and the multi-GeV p-\ike are shown separately for FC and PC events. The hatched region 
shows the Monte Carlo expectation for no oscillations normalized to the data live-time with statistical errors. 
The bold line is the best-fit expectation for f M <-> v T oscillations with the overall flux normalization fitted as a 
free parameter. 
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FIG. 4. The ratio of the number of FC data 
events to FC Monte Carlo events versus reconstructed 
L/E„. The points show the ratio of observed data to 
MC expectation in the absence of oscillations. The 
dashed lines show the expected shape for <-> v T at 
Am 2 = 2.2 x l(T 3 eV 2 and sin 2 29 = 1. The slight 
L/E v dependence for e-like events is due to contami- 
nation (2-7%) of CC interactions. 



n i i i i in i i i i 1 1 1 1 r 



1.5 



O 



c 
o 



D 



0.5 - 



• e-like 
o p.-like 

J I I I L 



9 



